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ABSTRACT 

We describe a web-based cgi calculator for constructing synthetic color- 
magnitude diagrams for a simple stellar population (SSP) using the Yonsei- Yale 
(YY) isochrone data base. This calculator is designed to be used interactively. 
It creates quick look CMD displays in (B-V) and (V-I) colors. Stochastic ef- 
fects on the CMDs are included. Output in tabular form is also provided for 
special purpose displays, or for combining the CMDs of different stellar popula- 
tions. This research tool has applications in studies of the stellar content of our 
Galaxy and external systems. It provides an easy way to interpret the CMDs 
in resolved stellar populations. It offers the means to explore the dependence 
of the integrated properties of unresolved stellar systems on stellar parameters 
(ages, chemical composition, binarity) and on the characteristics of their parent 
population (IMF slope and mass range). 



1. Introduction 

This program creates a synthetic color-magnitude diagram (CMD) for a simple stellar 
population (SSP), i.e. for a stellar population of a given age and chemical composition, 
that obeys a specified initial mass function (IMF) based on the prescription by Park & Lee 
(1997). Stochastic effects on the CMD are included. As an additional feature, the presence 
of a population of unresolved binaries can be taken into account (Woo et al. 2003). 

The code is designed for use with the Yonsei- Yale (YY) isochrones, and covers the pre- 
main sequence and hydrogen burning phases of evolution. The YY isochrones and luminosity 
functions are described in detail in papers by Yi et al. (2001, 2003; Paper 1 and Paper 3) 



-2- 



and Kim et al. (2002; Paper 2). Evolutionary tracks are available in Demarque et al. (2004; 
Paper 4). 

The YY isochrone tables cover metallicities Z (heavy element content by mass) in 
the range 0.0 to 0.08, and values of a enhancement corresponding to [a/Fe] = 0.0, 0.3 
and 0.6. Helium diffusion and convective core overshoot have been taken into account in 
calculating the evolutionary tracks. This first set of isochrones (Paper 1) was for the scaled 
solar mixture. The completing sets (for twice and four-times a-enhanced) were released in 
Paper 2. Two significant features of these isochrones are that (1) the stellar models start their 
evolution from the pre-main sequence birthline instead of from the zero-age main sequence, 
and (2) the color transformation has been performed using both the tables of Lejeune et al. 
(1998), and the older, but now modified, Green et al. tables (1987). 

The web interface to the CMD calculator can be found at the two YY web sites: 

http: / /www.astro. yale.edu/demarque/yyiso. html 

or 

http:/ /csaweb.yonsei.ac.kr/ kim/yyiso.html 



In addition, these two web sites provide links to the four YY papers listed above, and 
codes to interpolate between isochrones in age and chemical composition within the available 
parameter range. 

2. CMD calculator input parameters 

A number of input parameters must be defined in the construction of synthetic CMDs. 
The list of input parameters on the web interface is given below. The output CMD can be 
retrieved either in tabular or graphical form. 

• Alpha 

This is the composition parameter [a/Fe], which denotes the relative abundance by 
number with respect to the Sun of a enriched nuclei in the chemical composition 
mixture. The notation is the usual logarithmic relative abundance with respect to the 
Sun. In this notation, the metallicity of a star is defined by 

[X/H] = log(N x /N H )+ - log(N x /N H ) Q (1) 

where Nx and Nh are the relative abundances by numbers of element X and of hy- 
drogen, respectively. The abundance ratio [a/Fe] is defined correspondingly. Input 
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a The scaled-solar abundance ratios of metals are taken 
from Grevesse & Noels 1993. 



Note. — The abundance of the elements in logarithmic 
scale, logN e i/NH + 12, where N e i is the abundance by 
number. 



-4- 



values of a must be in the range 0.0 to 0.6 (by definition, a = corresponds to the 
solar mixture). 

• Metallicity 

The metallicity parameter Z for the synthetic CMD. Z is the mass fraction of the 
heavy elements as usually defined in stellar interior calculations (i.e. X + Y+ Z = 
1, where X and Y are the mass fractions of H and He). The parameter Z should be 
in the range (near) 0.0 to 0.08 (Note that Z = is not allowed). Table 2 provides a 
conversion between Z and [Fe/H] for different values of [a/Fe\. 

• Age 

The age of the stellar population (in Gyr) for the calculated synthetic CMD. 

• Number of stars 

The total number of stars N to be used in the calculation of the synthetic CMD. A 
maximum of 100,000 stars is allowed for the value of N. 

The total number of stars N in the synthetic CMD is defined as: 



where x is the IMF slope, and M upper and M\ ower define the range of masses included 
in the integration. The definition of x is the same as in the review by Tinsley (1980). 
The YY luminosity functions are tabulated for x = -1, 1.35 and 3. The x = 1.35 power 
law corresponds to the Salpeter IMF (1955). C is a normalization constant. The YY 
tables are normalized to 1000 stars in the mass range 0.5-1.0 M Q . 

• IMF slope 

Chosen value of the IMF slope x defined in eq. (2) (Note that x^O). 

• IMF upper limit 

The upper mass limit M upper of the IMF in eq.(2), in solar units. 

• IMF lower limit 

The lower mass limit Mi ower of the IMF in eq.(2), in solar units. 



The random number generator seed to be used for the synthetic CMD calculations. 
Varying ISEED illustrates the importance of stochastic effects in CMDs for a given 
set of input parameters. A knowledge of ISEED allows the user to retrieve previously 
made synthetic CMDs. 




(2) 



• ISEED 
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• Binary fraction 

The fraction of binary stars for the calculated synthetic CMD, given as a number 
between 0.0 and 1.0. 

• qmass 

Mass fraction for binary stars, given as a number between 0.0 and 1.0. 

3. Examples 

A series of examples are now presented that illustrate the effect of varying the input 
parameters one at a time. Fig. 1 shows the reference synthetic CMD against which other 
CMD examples can be compared, The following input reference parameters were adopted: 



REFERENCE PARAMETERS 
alpha - 0. 
Z = 0.004 
Age =8.0 
N = 1000 
IMF slope = 1.35 
-M U pp er — 5 . 
Mi ower = 0.6 
ISEED =1807 
Binary fraction = 0. 
qmass = 0.5 

Note that in every case below, only the one parameter discussed was varied; all other pa- 
rameters are as given in the above reference parameter list. 

• Varying chemical composition 

Metallicity Z 

Fig. 2 shows the effect of changing the metallicity Z from 0.004 to 0.02, all other input 
parameters being kept the same as in Fig. 1 . The reference heavy element solar mixture 
is the solar mixture of Grevesse & Noels (1993). It is given in Table 1, taken from 
Paper 2. 
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The parameter alpha 
Fig. 3 shows the effect of varying [a/Fe] from 0. to 0.3. Table 2, taken from Paper 2, 
lists the corresponding changes in total Z. The reader is referred to Paper 2 for a 
discussion of the different effects of variations in [a / Fe] in low-Z and high-Z mixtures. 

Relative H and He abundances 
Note that in the YY isochrone tables, the helium content by mass Y is kept fixed for 
a given value of Z. In all examples given here, (AY/AZ) = 2 has been assumed. This 
corresponds to Y = 0.238 for Z = 0.004, and to Y = 0.27 for Z = 0.02. 

• Varying the age 

In Fig. 4, the age has been changed from 8.0 to 5.0 Gyr. 

• Varying the total number of stars N 

In Fig. 5, the value of N is changed from 1000 to 500. As expected, stochastic fluctu- 
ations become larger as N decreases. 

• Varying the mass function parameters 

The distribution of stellar masses in the synthetic model depends on three parameters, 
the IMF slope x and the upper and lower masses, M upper and M [ower respectively. 
The effect of changing the IMF slope from x = 1.35 to 0.5 is shown in Fig. 6. The 
effects of changing independently the upper and lower mass limits to 3.OM and O.7M 
are shown in Fig. 7 and Fig. 8, respectively, all other parameters are as given in the 
reference parameter list. 

• ISEED 

Varying ISEED highlights the role of stochastic effects. For small N, stochastic effects 
become significant in rapid evolutionary phases. For example, the two CMDs shown 
in Fig. 9 and Fig. 10 are based on the same set of input parameters (except that 
ISEED=299 and 128, respectively). This introduces an uncertainty in the integrated 
energy distribution of the stellar population. 

• Varying the binary star content 

The presence of unresolved binaries can also affect the CMD significantly. Fig. 1 shows 
a synthetic CMD that does not include binaries. Fig. 11 and Fig. 12 are plotted 
for a binary fraction of 0.2, but with the different mass fractions qmass = 0.3 and 
0.7, respectively. The synthetic cluster CMDs by Woo et al. (2000) illustrate the 
importance of including binary systems, particularly in moderately rich stellar systems. 
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4. Conclusions 

This note describes a web-based CMD calculator designed for the study of stellar popu- 
lations in the Galaxy and distant stellar systems. This web based calculator, which is based 
on the YY isochrones, is a useful tool for studies of resolved star clusters, and for explo- 
rations of the sensitivity of the integrated light of stellar systems to stellar and population 
parameters. 

Future versions of this calculator will be based on a more comprehensive YY database 
covering an extended range of helium abundances, the inclusion of helium burning phases 
of evolution, often very significant in the CMDs of stellar systems, and the availability of a 
simple spectral energy distribution (SED) library. 

This research was supported in part by NASA grant HST-GO-10505.03-A. 
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Table 2. Conversion from [Fe/H] to Z 



[Fe/H] 


[a/Fe]=0.0 


Z 

[a/Fe]=+0.3 


[a/Fe]=+0.6 


-3.0 


0.000019 


0.000032 


0.000058 
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0.000062 
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0.000195 


0.000321 
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-1.5 


0.000615 


0.001012 


0.001807 


-1.0 


0.001935 


0.003174 


0.005627 


-0.5 


0.006021 


0.009774 


0.016990 


0.0 


0.018120 


0.028557 


0.047000 


0.5 


0.049711 


0.072793 


0.106471 


1.0 


0.110798 


0.142689 


0.177489 



